Anisotropic Superexchange for nearest and next nearest coppers 
in chain, ladder and lamellar cuprates 
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We present a detailed calculation of the magnetic couplings between nearest-neighbor and next- 
nearest-neighbor coppers in the edge-sharing geometry, ubiquitous in many cuprates. In this ge- 
ometry, the interaction between nearest neighbor coppers is mediated via two oxygens, and the 
Cu-O-Cu angle is close to 90°. The derivation is based on a perturbation expansion of a general 
Hubbard Hamiltonian, and produces numerical estimates for the various magnetic energies. In par- 
ticular we find the dependence of the anisotropy energies on the angular deviation away from the 90° 
geometry of the Cu-O-Cu bonds. Our results are required for the correct analysis of the magnetic 
structure of various chain, ladder and lamellar cuprates. 
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I. INTRODUCTION 

The magnetic interactions in the copper oxides are be- 
lieved to be governed by kinetic superexchange through 
the intervening oxygens. In tetragonal symmetry, one 
may view the CuO planes as consisting of clusters of 
four oxygens, forming a square whose center is occu- 
pied by a copper ion. These squares can be lined up 
along their edges, and then the nearest neighbor (NN) 
Cu-O-Cu bond makes an almost 90° angle. Another 
ubiquitous configuration is formed when the squares are 
connected along their corners, in which case the near- 
est neighbor Cu-O-Cu bond is linear, having an angle of 
180°. Typical examples for edge-sharing compounds are 
La6CasCu2404il!J (where the. angle is 91°) and CuGeOa 
(where the angle is s=s 98° )E Corner-sharing configura- 
tions characterize the copper oxide planes in the parent 
compounds of the higrtjT c cuprates, and the chains in 
Sr 2 CuO30 and SrCu02l3 Some compoundsJnclude both 
types of bonds, for example Si^CusO-jC^.B The various 
Cu-Cu bond geometries in this material are the same 
as those for the nearest-neighbors and nextppearest- 
neighbors (NNN) in the chains in Sri4Cu2404iB and the 
interladder bonds in Sr„_iCu„ + i02 n -Q 

The magnitude and the sign of the magnetic inter- 
actions in the two types of bonds are expected to be 
quite different. According to the set-called Goodenough- 
Kanamori-Anderson (GKA) rulesjj the leading isotropic 
superexchange of a 180° bond between two magnetic ions 
with partially filled d shells is strongly antiferromagnetic, 
while the leading order of a 90° superexchange is ferro- 
magnetic, and much weaker. In the Cu-0 case, the rea- 
son for this is that for the corner-sharing geometry, the 
2p(j orbital hybridizes with the two neighboring Cu ions, 
yielding a significant contribution to the kinetic superex- 
change (which is antiferromagnetic). In contrast, in the 
edge-sharing configuration the 2p CT orbital on the oxy- 
gen, which hybridizes with a 3d orbital on one copper, 
is almost orthogonal to that 3d orbital on the nearest- 
neighbor Cu ion, thus blocking the antiferromagnetic su- 
perexchange via a single oxygen. The leading magnetic 



coupling in this case is given by the next order perturba- 
tion terms, and is therefore weaker and of the opposite 
sign. 

Higher-order perturbation terms also determine the 
magnetic anisotropics, in both types of bonds. These 
anisotropics are responsible for various observable quan- 
tities, like the gaps in the spin wave spectrum, the spin 
orientations in space, etc, and hence are of much inter- 
est. The magnetic couplings of the linear Cu-O-Cu bond 
were investigated in great detail (see Refs. g and [10| and 
references therein), yielding the in-plane and the out-of- 
plane gaps of the family of compounds with structures 
similar to that of La2Cu04, as well as the antisymmetric 
Dzyaloshinskii-Moriya interaction in the orthorhombic 
phase. In particular, the role of the on-site Coulomb ex- 
change and the spin-orbit interaction in producing the 
various anisotropies was clarified in tetragonal and or- 
thorhombic symmetries. 

In the almost-90° bond (see Fig. [I]) the leading order 
magnetic exchange is small. Therefore, higher-order per- 
turbation processes, as welL-as-details of the structure like 
the presence of side groupsllilll3 have a significant contri- 
bution. For the same reason, next-nearest-neighbor Cu- 
Cu couplings are expected to be much more important 
than in the case of the linear Cu-Cu bond configuration. 

Previous discussions of this geometry include an anal- 
ysis of the dependence of the leading nearest-neighbor 
isotropic coupling on the small angular deviation away 
from 90°, (5,Ej££l which has been found to be domi- 
nated by the_pn-site Coulomb exchange interaction on 
the oxygens,Ej and the nonlocal exchange between the 
coppers and the oxygens .E31L3 The magnetic anisotropies 
have been calculated only for the nearest-neighbor, 
strictly 90°-bond, by Yushankhai and Hayn.ll3 We com- 
pare below their results with ours. 

The aim of this paper is to present a detailed calcula- 
tion of both the isotropic and the anisotropic magnetic 
interactions in the nearly-90° configuration shown in Fig. 
Fy, for nearest-neighbor and next-nearest-neighbor cop- 
per ions. Our calculation is based on the perturbation ex- 
pansion of a Hubbard model around the half- filled ground 



state, in which there is a single 3d hole on each copper 
ion, and the oxygen 2p states are completely full. Since 
the spin of the Cu hole is arbitrary, this ground state 
is 2 N -fold degenerate, where N is the number of cop- 
per ions. The superexchange magnetic Hamiltonian is 
obtained as the effective interaction within this degen- 
erate manifold. The microscopic Hamiltonian we con- 
sider includes hopping between all orbitals on the cop- 
per d-states and on the oxygen p-states, and between 
the p states on neighboring oxygens. The Hamiltonian 
also contains the spin-orbit interactions on the copper. 
Those on the oxygen are much weaker and are therefore 
neglected. We also include all local Coulomb interac- 
tions on the copper and on the oxygen, and the non-local 
Coulomb exchange between the copper and the oxygen. 
General expressions for the effective magnetic Hamil- 
tonian, based oh. this micorscopic Hamiltonian have been 
derived before£rli2l In order to keep the paper self- 
contained we reproduce in Sec. II the main steps of the 
derivation. We find in Sec. II that the magnetic Hamil- 
tonian has the form 



n 



(ij) l 



M cMcM 
NN D i D j 



tM CM C/A 

J NNN a i D j )' 



(1) 



in which \i denotes the Cartesian component of the spin 
and Jnn and Jnnn are the magnetic couplings between 
nearest-neighbors and next-nearest-neighbors, respec- 
tively. It is convenient to define the coordinate system 
for the spin components such that the x and y directions 
are in the Cu-0 plane, along the bonds between the cop- 
pers, and the z direction is perpendicular to the plane. 
The leading magnetic coupling for both NN and NNN is 



r 



J x + J y + J z 



(2) 



The anisotropic couplings are then naturally given by 
J op , for the out-of-plane anisotropy, and J pd , for the 
in-plane one: 



Top _ jz _ J "b •[_ rpd _ J J 



(3) 



The notation "pd" stands for pseudo-dipolar, see Ref. 



The parameters that determine the magnitude and the 
sign of the magnetic couplings are the Cu-0 and 0-0 
hopping matrix elements, the on-site (single particle) en- 
ergies on the oxygen and on the copper, the spin-orbit 
coupling constant A and the various Coulomb matrix el- 
ements. The latter are parametrizedO in terms of the 
Racah parameters into the on-site leading order interac- 
tions, and the residual remaining interactions, which are 
small. However, they, as well as A, are necessary for the 
generation of the magnetic anisotropies. All these param- 
eters depend on the crystal symmetry, and hence on the 
angle 8 (see Fig. [I]). Adopting the plausible assumption 
that the most important sensitivity to small deviations 



from 90° occurs in the Cu-0 hopping matrix elements, 
we have included only their dependence on the angleJiS 
and calculated the angular dependence of the magnetic 
interactions for small angles 8. The explicit expressions 
for the magnetic couplings are given in Sees. Ill and IV. 
Here we summarize the results, which are depicted in Fig. 
g for the nearest-neighbor couplings, and in Fig. g for 
the next-nearest-neighbor ones. 

The numerical estimates are computed using the fol- 
lowing parameters. We take the on-site energies on the 
oxygen to be e Px = e p =3eV and e Pi =2eV.E!J Those 
on the copper are assumed for simplicity to be identk 
cal, e a =1.8eV (a = 0, l,x, y, Z, see definitions below) .0 
The spin-orbit coupling on the copper is taken to be 
A=0.1eV. The on-site Coulomb matrix elements ne- 
cessitate the Racah parameters A, B, and C for the 
copper, and Fq and F2 on the oxygen.0tL9 These are 
chosen as A=7.0eV, B=0.15eV, C=0.58eV, F =3.1eV, 
and F2=0.28eV. There is no reliable estimate for the 
non-local Coulom b exchange between the copper and 
the oxygen EJO'til We therefore take as a representa- 
tive estimate a value in the range between 0.02eV and 

o.iev.yyy 

The hopping matrix elements can be .expressed 
in terms of the Slater-Koster parameters JLM to = 
—V3(pda)/2, t\ = —(pda)/2, and i 2 = (pdir), for 
the Cu-0 ones, and £3 = (l/2)((ppcr) + (ppn)), £4 = 
(l/2)((ppa) — (ppir)), and £5 = (ppir) for the 0-0 matrix 
elements. We have used the values (pda) = 1.5eVLrE2l 
and (ppn) = — 0.6eVll3, and used the relations (ppn) — 
- \ {ppcr) and {pdn) = - \ {pdaJB. 

Figure || depicts the angular dependence of J^ N , ^NN' 
and J^n- Th e three curves in each figure are obtained 
by choosing different representative values for the non- 
local Cu-0 Coulomb exchange matrix element, K. At 
strictly 90°, J^n i s negative (ferromagnetic) and small, 
w -0.02eV (-0.04eV, -0.07eV) for K = 0.02eV (0.05eV, 
O.leV). Its value is determined mainly by the resid- 
ual Coulomb interactions. As the angle deviates from 
90°, J§n approaches zero and changes its sign at about 
8 w0.05 (0.065, 0.09) (2.9°, 3.7°, 5.2°). These results 
agree with those found before in Refs. [H 



and 
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The 



out-of-plane anisotropy J^n ^ s relatively large, and neg- 
ative, in agreement with the findings of Ref. p6|. Like all 
other anisotropies, its magnitude is proportional to A 2 . 
However, for delicate reasons related to "ring-exchange" 
processes (see below), the Coulomb matrix element that 
scales its magnitude is the on-site interaction on the oxy- 
gen, leading to its comparatively high value, ~ — 1.3meV 
at 90°. (The out-of-plane anisotropy increases slightly 
with increasing K.) The nearest-neighbor in-plane 
anisotropy J|J N is scaled by the small residual Coulomb 
interactions. It vanishes at 8 — 0, and stays quite small 
away from that value, varying approximately linearly 
with 8, « (-0.1 to -1.7 K)8 meV. 

The case of an ideal 90° nearest-neighbor bond has 
been recently discussed in Ref. [ll| These authors have 



specialized to materials of the type A2CU3O4CI2, with 
A=Ba or Sr. They have neglected the non-constant on- 
site Coulomb interactions on the oxygen and the non- 
local Cu-0 Coulomb interaction, but have taken into ac- 
count the local orthorhombic symmetry, by allowing the 
Cu on-site energies e x and e y to be different. They there- 
fore obtained a small in-plane anisotropy, of the order of 
0.2^eV. 

The analogous results for the next-nearest-neighbor 
couplings are summarized in Fig. H. These necessitate 
additional perturbation processes, which involve hopping 
between nearest-neighbor oxygens. We find that J^nn 
is about 20meV at 90°, and has a smooth linear depen- 
dence on 5 away from it, remaining antiferromagnetic for 
small angles, in agreement with the findings of Refs. n3 
and |l|. 

As in the case of the nearest neighbors, also the 
anisotropic coupling J^nn * s relatively large and nega- 
tive, being rs — 0.036meV at 90°, while J^nn IS extremely 
minute, i=a — 6/xeV. The out-of-plane anisotropy is again 
dominated by the "ring-exchange" processes and the in- 
plane anisotropy by the small residual Coulomb interac- 
tions AU. 

Obviously, the results summarized in Figs. and 
depend on the details of the parameters, e.g., the hop- 
ping matrix elements or the Coulomb Racah coefficients. 
The remaining Sections of the paper are devoted to a de- 
tailed discussion of the derivation and the choice of these 
parameters. 



orbital angular momentum vector between the crystal- 
field states a and /3. The non-zero matrix elements 



are Ln 



-2», LI 



L 0y 



T X 



~ L ly 



y/Ei, 



L v zx = ~L x zy — L xy — i, and L* a „ — Lp a . The last term 
in (||) is the Coulomb interaction, with U a p^s = {o-5\[3"f) 
in the notations of Table A26 in Ref. [17| The Hamilto- 
nian of the oxygen ions is 



"Ho = ^2 £ n p\ 



7 j KPnmPqnaPq'ma + " c j 



qna 



qq' 



r, 2-~i ^ n i n 2n3n 4 P qnia P qn2(7 iPqn 3 a'Pqn i (7, (6) 



in which p qna creates a hole in one of the three p or- 
bitals, p x , p y , and p z (denoted by n) on the oxygen at 
site q, with energy e„. The second term in (||) describes 
the hopping between the O-ions, and the last term is 
the Coulomb interaction on the O-ions. Finally, TLcu-o 
describes the kinetic energy of hopping between the Cu 
and the O ions, and the Coulomb exchange interaction 
between them, 

Hcu-O = Z~2{tn a P\nA^a + he) 

iqa 

an 
+ Yl K <*mf3n4 a(7 pl ma , diffa' Pqna- (?) 

a(3mn 



II. THE MAGNETIC HAMILTONIAN 

As is discussed above, the magnetic Hamiltonian is de- 
rived from a microscopic Hamiltonian. The latter can be 
written as folllows: 



H = Hcu + 'Ho + Wc«- 



O) 



(4) 



A significant simplification of the perturbation expan- 
sion is achieved by first treating the spin-orbit interac- 
tions exactly, leaving the expansion-in orders of the spin- 
orbit coupling, A, to the final stage.Q This is accomplished 
by introducing the unitary transformation which diago- 
nalizes the single-particle part of Ticu 



diaa — / X rn aa\aa' C iaa'i 



(8) 



with obvious notations. Explicitly, the Cu-ion Hamilto- 
nian is 

H C u = 2_, e ad iaa d ia(T + - 22 L «/3 • [e\<j<j'd iacr d l(itjl 



+ „ 2-*/ UaPjSUiaudj^p^d^cridi&o 



(5) 



a{3~y8 



where d\ a(J creates a hole with spin a in the crystal-field 
state a at site i, of site energy e a . For tetragonal sym- 
metry we label the crystal-field states as |0) ~ x 2 — y 2 , 
|1) ~ 3z 2 — r 2 , \z) ~ xy, \x) ~ yz, and \y) ~ zx, where the 
z axis is perpendicular to the plane and |0) is the lowest 
energy single-particle state. The second term in (||) is the 
spin-orbit interaction, where A is the spin-orbit coupling 
constant and \j a p denotes the matrix elements of the 



where c] ocr creates a hole in the exact eigenstate a of the 
Hamiltonian which consists of the crystal-field and the 
spin-orbit interaction on the copper. These states have 
a site label i, a state label a, and a pseudo-spin index a. 
One then has 



A 



7 , e ad iacr d laa + — 2_^ L Q/3 • [o]a<r'd iarr d. 

iaa iap 

aa' 

— 2_^, ^aC iaa Ci aai 



where E a and \m ao \ aa i are determined by 

Eb[m,j b]crcr' — £j[m*fb]cTCT' 
0a t 



i/3a' 



(9) 



(10) 



witl1 E Q m /3a( m <m) t W = S a p$aa>- When A -> 0, each 
state | a) approaches one of the states |a). Using this def- 
inition, the index a runs over the values 0, 1, z, x, and 
y. A detailed discussion of this transformation is given 
in Rcf. | 

To apply the pertubation expansion, we divide the 
Hamiltonian TL into an unperturbed part, Ho, and a per- 
turbation term TL\. The part 7io contains the single 
particle Hamiltonians on the coppers and on the oxygens, 
and the leading on-site Coulomb potentials. The pertur- 
bation Hamiltonian contains the kinetic energy and the 
residual Coulomb interactions. As is known,£3 the on- 
site Coulomb interactions can be parametrized in terms 
of the Racah coefficients. In tetragonal site symmetry, 
those on the copper are parametrized by the Racah pa- 
rameters A, B, and C, with A ^S> B and A 3> C, 
and those on the oxygen by Fb and F%, with Fq 3> F^. 
We choose the on-site leading Coulomb interactions to 
be U = U aaaa = A + AB + 3C for the copper, and 
U q = Unnnn = Fq + 4F 2 on the oxygen. Consequently 
the unperturbed Hamiltonian is 



with 



U 



^0 — / v ■&aCia (T Cia<T + r, / J c iacr C iba' Ci ' 



bo' Cii 



u, 



/ , tnPqngPqno + „ / jPgnaPgn' cs'P<l nl a 'P<t na - (llj 



qua 



The perturbation Hamiltonian is 

Wi = Tihop + AHc, 
in which the hopping term is 

Hhop — / \ \kan\a' tjCj a<y iPqnt7 + MC) 

iqan 

era' 
> / A nmPqno-Pq' mo- + " c j- 
qq> 



(12) 



(13) 



Because of the transformation rf|), the Cu-0 hopping 
becomes spin-dependent 



[tanlcr'a — / J * an [Tn aa \ 



(14) 



The term AHc contains the (small) additional on-site 
Coulomb interactions, and the non-local Cu-0 Coulomb 
potential, 

A?ic = o E AU ^28 3 s i (abcd)4 asi 4 bs2 c ic83 c idSi 



iabcd 
S1S2S3S4 



cy / 4 ^Un 1 n 2 nzn 4: 'Pqn 1 a'Pqn2<T'P ( l n 3v' 



Pqn^a 



nin2n^ri4 



^2 K sa' 8 'o{ambn)c\ as p < (ima ,Ci bal p qna , (15) 



iqabmn 

ss' aa' 



AU SlS2S3Si (abcd) 

Af/ Q/ 3 7 «5[(m Qa ) t m,5 d ] SlS4 [(m /3h ) t m 



'S1S2S3S4 

E 



and 



K sa , s t a {ambn) 



2_, Kam,f3n[m a a] as [m{3b]cr' s' 
a/3 



(16) 



(17) 



Here we have defined AU, 



fc^t-/ nn' n' n ■ ^n 



t/„ 



a/3 /3a 



n\ ^ rii or or 712 7^ n.3, AC/ 



= ^ q /3/3q - U Q and 
For a ^ or /? ^ 7, and 
C/ripvolves only the small 
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Racah coefficients _B, C, and F2O 

All the perturbation contributions resulting from Hi 
begin and end within the 2 N -fold degenerate ground- 
state manifold of Ho, each state of which has one hole at 
each copper site, with arbitrary spin a. We will denote by 
"0" the ground state of the single-particle Hamiltonian 
on the copper, and will take its site energy to be zero. 
For the sake of clarity, we divide the perturbation chan- 
nels into three groups. Group "a" includes the processes 
in which there are two holes on the copper in the inter- 
mediate step; group "b" includes those in which there 
are two holes on the oxygen in the intermediate step, or 
processes where the two holes exchange their correspond- 
ing coppers by going around the ring formed by the two 
coppers and the intervening oxygen (for example, cop- 
pers i and j, and oxygens q and q' in Fig. |I|); group "c" 
contains the contribution from the non-local Coulomb 
exchange between the coppers and the oxygens B3 

The contribution of channel "a" to the magnetic cou- 
pling of coppers i and j is 



E 

ab 

otft-yS 



(U + E a ){U + E b ) 



(Tr{a ■ S t f ^(m ab y mSa f^}Tr{<j ■ S J (m /30 )tm 70 } 

-Trjcr • S i T ^(m Q!6 ) t m 5 ocr • S i (m /30 ) t m 7a 'f a : 'j|j 

+ (* - j), (18) 

in which Sj is the spin on the copper at site i in the new 
orbital ground state \a) = |0), 



Si — 77 7 , c iOffMjff ,c iO<r'] 



(19) 



and the "Tr" are carried out in spin space. E3 We have in- 
troduced in ( jig ) the notation 7"^ for the effective matrix 
element for hopping from state a on copper j to state b 
on copper i. These are different in the case where the 
two coppers are nearest neighbors, and when they are 



next-nearest neighbors. In the first case such a process 
can be achieved through a single oxygen, yielding 



[■*ba\ c ''<7 ~ y. , rbnlv'virna 
qnai 



(20) 



to lowest possible order in perturbation theory. Figure 
^ depicts the direct hopping from copper i to oxygen q, 
together with the two possible indirect hoppings, going 
through oxygen q' and oxygen q" . To account for these 
processes one has to use the following replacement in (pp|) 



\PnaW\t 



rnajcirr 



/ t t nm lt ma \a lC r- (21J 

inq 



When the two coppers are NNN, the 0-0 hopping is 
essential for bringing the two holes to the same copper. 
(For example, the bond j'-j in Fig. requires the q-q" 
hopping.) In this case we have 



where 



i^ba\ cr ' c ' 



[ft] 



V- 1 - 
e 



/ , c \Pbn\ a ' a i\- na\oioi 
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E 

mq' 



.+11 \+1 i 1 



(22) 



(23) 



The perturbation contributions coming from channel 
b, when coppers i and j are NN, yield 



Hb(i,j) 



6nn>- 



ranqq' 
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u„ 



+ (i-V)- 



Tr 



[ . alii' fl'i . c .pi Zqi \ 



E 



At/, 



mm nn 



{U q + e n + e n ')(U q + e m + e m <) 
(Tr{a ■ S^J^}Tr{a ■ S 3 -t* ,*&} 



(24) 



The term with the (1 — 5 qq i) in front arises in the chain 
geometry and does not have-tie on-site Coulomb inter- 
action in the denominator.EJEil As non-local Coulomb 
interactions between the oxygens are ignored here, there 
is no analogous contribution in the second sum of (24). 
When the two coppers are NNN, one must invoke the 
0-0 hopping. We then find that each pair of matrix 
elements t has to be replaced as follows 



+1^+13 



rpqifqj _ fiirpqj 



(25) 



where both i and T are matrices in spin space, given by 
Eqs. (|14| ) and (23), respectively. 
Finally, channel c gives 



n c (i,j) 



v a.m-yn 






tn^m 

+(* <-► j), 



Tr[a ■ SiCmaoW* • S^m^o} 



(26) 



when the two coppers are NN. The corresponding expres- 
sion for NNN coppers is obtained from (pi) by replacing 
each matrix element t by T. 

The total effective magnetic interaction is the sum of 
the three groups. This is now expanded up to second 
order in the spin-orbit coupling A. This requires the ma- 
trix to, Eq. (ph, up to first order in A only. Using Eq. 
(|lO|) we have 



m 



^Qfl"0 



2 t 



E a 



(27) 



A further significant simplification of the expressions 
is achieved when one takes into account the_symmctry 
properties of the Coulomb matrix elementsmLB AU a pj5 
vanishes unless the values of a, /3, 7, and 5 are such 
that the products <r(a)a(5) and cr(f3)a(j) are propor- 
tional to each other. Here we use the convention that 
a(a = 0) and a(a = 1) are the unit matrix. Similarly, 
the non-vanishing matrix elements of At/ ni „ 2n3 „ 4 satisfy 
a(ni)a(ri4) ex a(n,2)<j{ns), and those of K am p n vanish 
unless a(a)a(n) ex er(/3)er(m). 

In the following we present the general expressions in 
the form 



lie, .;n ( J - ^Trr J s, • s s + s,rs 



3' 



(28) 



where J includes all contributions to zeroth order in 
A, while the matrix T contains the contributions which 
necessitate the spin-orbit interaction, and is therefore 
second-order in A. The application of these general ex- 
pressions to the specific Cu-Cu bonds will be carried out 
in the next Section. 

(a) Channel a [Eq. (|§)] yields 
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(b) Channel b [Eq. (|§)] yields 
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and 
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(c) Channel c [Eq. @] yields 
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(36) 



For simplicity, we have written the results for J and 
r, Eqs. (J2E-36) for the nearest-neighbor Cu-Cu bond. 
The analogous expressions for the next-NN bond are ob- 
tained using the replacements (E3) and (25). 



III. THE MAGNETIC COUPLINGS 



A. The nearest— neighbor 90° bond 



We list in Tables S and pi the hopping matrix elements 
between Cu and O for the 90° configuration (see the In- 
troduction and Fig. f| for the notations). Using these 
values, we find that the leading order contributions to 
the magnetic couplings come from channels "b" and "c" , 
with 



J — Jb + Jc — 

167J A[/p x p y p xP!; ^ 8/J 
el(2e p «+U q y £**.<*.■ 



(37) 



in accordance with the GKA rules. Inserting 

the numerical values of the parameters, we find 
J sa — 8meV— OJTXop^Op^ The leading order magnetic 
anisotropy in this case is the out-of-plane one, 



leading order 
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(38) 



with T zz w — 1.3meV. The processes yielding the latter 
are depicted in Fig. [| 

The remaining small anisotropics resulting from the 
on-site Coulomb potential on the oxygen, and from the 
non-local Coulomb exchange between the copper and the 
oxygen, are listed below. We express those in the coordi- 
nate system depicted in Fig. ||. To obtain the couplings 
in the coordinate system discussed in the Introduction, 
one has to rotate by 45°. 
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(39) 



To obtain these results we have used the relationsEZI 
Caoio = 4B + C, UxQxo = U y o y o — 3B + C, 
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K 0py op y = -Kq Px q 1l .,K ZPxZPx = K zPyZPy . In addition, 



from Tables || and O the only non-zero effective Cu-Cu 
hoppings are 



j-U _ M _ 
l z\ ~ l z\ — 
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(40) 



Numerical estimates of these expressions yield Ar| z rs 
40//eV, Ar^ w 7/xeV, Arg^ « -9/xeV, Arf « -if x 
0.003, Arjf w -A" x 0.001. 



B. The next— nearest— neighbor 90° bond 
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Here T 2 ^ = r j,x . As before, this coordinate system has to 
be rotated by 45° in order to produce the magnetic cou- 
plings in the form discussed in the Introduction. That 

is, 
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(46) 



Tables III and IV list the hopping matrix elements 
T lq n for a Cu-O-O-Cu process in the 90° conhgurtion 
(the notations are shown in Fig. 0). In these Tables, 

*3 = tpxPx = '■'PyPy ' 4 = tp x p y = _ *Pa,p« ) an Q *5 = tpzPz ■ 

The effective hopping matrix elements between two cop- 
pers, 
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which do not vanish arc 
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(42) 



As opposed to the NN-bond, in this case hopping be- 
tween the ground state orbitals of NNN coppers is pos- 
sible, via the two oxygen orbitals p x and p y , which are 
connected by £4, see Fig. [| We find contributions to the 
coupling J from all three channels, 

J = J a + Jb + J ci 



Ja = 



4* U ° 



Jb 
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with J i=s 0.02eV The leading order anisotropy is the 
out-of-plane one, and comes mainly from channel "b" 
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leading order 
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(44) 



r kfadin g order ~ -30^eV. The remaining non-diagonal 
anisotropies, caluclated in the coordinate system of Fig. 
p, are 



IV. SMALL DEVIATION FROM 90° 

To calculate the magnetic couplings as function of the 
angle S (see Figs, [l] and 0) we use the following forms for 
the hopping matrix elements 

t 0px — t cos (5(1 — 2 cos 2 5) + 2*2 sin 5 cos S, 

t^ = ~*o sin 5(1 - 2 cos 2 6) + 2* 2 sin S cos 2 S, (47) 
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with *?„ = *?„ and t q r , 
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and 



4L = ~*2 sin (5, 

tl q Pz =t 2 cos6, (50) 



with £«i = tg 2 



For the next-nearest neighbor hopping we need the 
following matrix elements: 



-Lp x a 1^3 *p x a ~r ^4^p y a) ; 

J-PyCr. ^S^PyQ "T t4tp x a J? 



for a = 0, 1,#, and 



-*-Pz a — h b v p z ct 
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(51) 



(52) 



for a = x,y. The resulting explicit expressions for the 
magnetic couplings are very long, and we therefore skip 
them. Instead, we have used the results above to pro- 
duce the curves in Figs. M and S, to obtain the magnetic 
couplings as function of the angular deviation S. 



V. DISCUSSION 



We have presented a detailed calculation of the mag- 
netic interaction between nearest-neighbor and next- 
nearest-neighbor coppers in the edge-sharing geometry, 
and obtained numerical estimates for the various cou- 
plings as function of the angular deviation from 90°. 
These numerical estimates are crucial for the analysis of 
the magnetic structures of many chain, ladder and lamel- 
lar cuprates. Our calculation is based on a perturbation 
expansion of a general Hubbard Hamiltonian. It has been 
found before that for the magnetic anisotropics of the lin- 
ear Cu-O-Cu bond, this expansion is quite reliable^ On 
the other hand, it has been argued, (again for the ISC- 
bond), that perturbation theory fails to yield reasonable 
values for the leading magnetic isotropic intractionsJEn 
because the hopping matrix elements t are not necessar- 
ily small compared with the on-site energies. However, 
the almost 90° case discussed here is different, because 
of the appearance of the small Coulomb matrix elements 
in the expansion. It therefore can be expected that the 
perturbation expansion for the present case yields reli- 
able estimates. Indeed, the comparison of our results for 
the NN isotropic, energy with those obtained from exact 
diagonalizationlU seem to support this conclusion. 

Our results show that the out-of-plane anisotropy is 
negative, both for NN and for NNN coppers. This in- 
dicates an easy axis perpendicular to the Cu-0 plane, 
in agreement with Ref. 111. We find that the pseudo- 
dipolar interaction between nearest-neighbors vanishes 
at strictly 90°, and is minute for a small deviation away 
from it. As has been shown in Ref. [l(| this result is 
modified when one allows for a difference between the 
Cu on-site energies e x and e y . It seems that this should 
be the case in materials like S^CusO/iC^, where some of 
the copper ions, lose their local tetragonal symmetry: It 
has been foundO that experimental data on Sr2Cu30.4Cl2 
imply a finite value for this energy. This means that 
the interpretation of the data necessitates the inclusion 
of such effects, or of dipolar interactions. In the same 
manner, it is expected that our numerical estimates as 
function of the angle S will be useful in the analysis of 
other cuprates. 
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FfG. f . Edge-sharing Cu-0 configuration. The angle of the Cu-O-Cu bond is n/2 — 28. Open circles denote oxygens and 
black circles are the Cu's. 
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FIG. 2. The nearest-neighbor magnetic couplings Jnn for K=0.02 eV (bold solid line), K=0.05 (solid line) and K=0.1 eV 
(dashed line). 
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FIG. 3. The next-nearest-neighbor magnetic couplings Jnnn for K=0.02 eV (bold solid line), K=0.05 (solid line) and 
K=0.1 eV (dashed line). 
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FIG. 4. The Cu-0 effective nearest-neighbor hopping in the 90° bond configuration, allowing for nearest-neighbor 0-0 
hopping. Open circles denote oxygens and black circles are the Cu's. The processes contributing to t' q are shown by arrows. 
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FIG. 5. Cu and orbitals which are involved in the processes leading to the out-of-plane anisotropy F 22 . Here and below, 
the shaded (white) area indicates positive (negative) phase. 
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FIG. 6. Next-nearest-neighbor hopping between two d x 2_ i orbitals via the p x and p y orbitals. 
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FIG. 7. Cu-O-Cu geometry for an angle ty/2 — 28. 



TABLE I. The Cu-O hopping matrix elements t lq n for the 9(f bond. Upper signs: t %q , lower signs: t 3q 
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TABLE II. The Cu-0 hopping matrix elements t„„ for the 90° bond. Upper signs: t 3q , lower signs: t l 
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TABLE III. The Cu-O-O hopping matrix elements T l a q n for the 90° bond. Upper signs: T iq , lower signs: T jq 
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TABLE IV. The Cu-O-0 hopping matrix elements T£ n for the 90° bond. Upper signs: T 3q , lower signs: T il 
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